Mon. Not. R. Astron. Soc. 000, mLuy (2000) Printed 1 February 2008 (MN WF&L style file vl.4) 



Gravitational waves from cosmological compact binaries 
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ABSTRACT 

We consider gravitational waves emitted by various populations of compact binaries 
at cosmological distances. We use population synthesis models to characterize the 
properties of double neutron stars, double black holes and double white dwarf binaries 
as well as white dwarf-neutron star, white dwarf-black hole and black hole-neutron 
star systems. 

We use the observationally determined cosmic star formation history to reconstruct 
the redshift distribution of these sources and their merging rate evolution. 
The gravitational signals emitted by each source during its early-inspiral phase add 
randomly to produce a stochastic background in the low frequency band with spectral 
strain amplitude between ~ 10 -18 Hz -1 ' 2 and ~ 5 x 10~ 17 Hz~ 1//2 at frequencies in 
the interval [~5x 1CT 6 - 5 x 1CT 5 ] Hz. 

The overall signal which, at frequencies above 10 Hz, is largely dominated by 
double white dwarf systems, might be detectable with LISA in the frequency range 
[1 — 10] mHz and acts like a confusion limited noise component which might limit the 
LISA sensitivity at frequencies above 1 mHz. 

Key words: gravitation - stars: formation - stars: binaries-gravitational waves. 
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1 INTRODUCTION 

Binaries with two compact stars are the most promising 
sources for gravitational radiation. The final phase of spi- 
ral in may be detected with ground-based (LIGO, VIRGO, 
GEO and TAMA) and space-borne laser interferometers 
(LISA). This has motivated researchers to model gravita- 
tional waveforms and to develop population synthesis codes 
to estimate the properties and formation rates of possible 
sources for gravitational wave radiation. 

Since there is not yet a single prescription for calculating 
the gravitational emission from a compact binary system, it 
is customary to divide the gravitational waveforms in two 
main pieces: the inspiral waveform, emitted before tidal dis- 
tortions become noticeable, and the coalescence waveform, 
emitted at higher frequencies during the epoch of distortion, 
tidal disruption and/or merger (Cutler et al. 1993). 

As the binary, driven by gravitational radiation reac- 
tion, spirals in, the frequency of the emitted wave increases 
until the last 3 orbital cycles prior to complete merger. 

With post-Newtonian expansions of the equations of 
motion for two point masses, the waveforms can be com- 



puted fairly accurately in the relatively long phase of spiral 
in (see, for a recent review, Rasio & Shapiro 2000 and refer- 
ences therein). Conversely, the gravitational waveform from 
the coalescence can only be obtained from extensive numer- 
ical calculations with a fully general relativistic treatment. 
Such calculations are now well underway (Brady, Creighton 
& Thorne 1998; Damour, Iyer & Sathyaprakash 1998; Rasio 
& Shapiro 1999). 

In this paper, we consider the low frequency signal from 
the early phase of the spiral in, which is of interest for space- 
borne antennas, such as LISA. For this purpose, we use the 
leading order expression for the single source emission spec- 
trum, obtained using the quadrupole approximation. Our 
analysis includes various populations of compact binary sys- 
tems: black hole-black hole (bh, bh), neutron star-neutron 
star (ns, ns), white dwarf-white dwarf (wd, wd) as well as 
mixed systems such as (ns, wd), (bh, wd) and (bh, ns). 

For some of these sources [(ns, ns), (wd, wd) and (ns, 
wd)], statistical information on the event rate can be inferred 
from electromagnetic observations. In particular, there are 
several observational estimates of the (ns, ns) merger rate 
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obtained from statistics of the known population of binary 
radio pulsars (Narayan, Piran & Shemi 1991; Phinney 1991). 

A rather large number of close white dwarf binaries have 
recently been found (see Maxted & Marsh 1999 and Moran 
1999). However, it is customary to constrain the (wd, wd) 
merger rate from the observed SNIa rate (see Postnov & 
Prokhorov 1998) . Also the population of binaries where a ra- 
dio pulsar is accompanied by a massive unseen white dwarf 
may be considerably higher than hitherto expected (Porte- 
gies Zwart & Yungelson 1999). 

Since most stars are members of binaries and the forma- 
tion rate of potential sources of gravitational waves may be 
abundant in the Galaxy, the gravitational-wave signal emit- 
ted by such binaries might produce a stochastic background. 
This possibility has been explored by various authors, start- 
ing from the earliest work of Mironovskij (1965) and Rosi 
& Zimmermann (1976) until the more recent investigations 
of Hils, Bender & Webbink (1990), Lipunov et al. (1995), 
Bender & Hils (1997), Giazotto, Bonazzola & Gourgoulhon 
(1997), Giampieri (1997), Postnov & Prokhorov (1998), and 
Nelemans, Portegies Zwart & Verbunt (1999). This back- 
ground, which acts like a noise component for the inter- 
ferometric detectors, has always been viewed as a potential 
obstacle for the detection of gravitational wave backgrounds 
coming from the early Universe. 

In this paper we extend the investigation of compact bi- 
nary systems to extragalactic distances, accounting for the 
binaries which have been formed since the onset of galaxy 
formation in the Universe. Following Ferrari, Matarrese & 
Schneider (1999a, 1999b: hereafter referred as FMSI and 
FMSII, respectively), we modulate the binary formation 
rate in the Universe with the cosmic star formation history 
derived from observations of field galaxies out to redshift 
z ~ 5 (see e.g. Madau, Pozzetti & Dickinson 1998b; Steidel 
et al. 1999). 

The magnitude and frequency distribution of the in- 
tegrated gravitational signal produced by the cosmological 
population of compact binaries is calculated from the dis- 
tribution of binary parameters (masses and types of both 
stars, orbital separations and eccentricities). These orbital 
parameters characterize the gravitational- wave signal which 
we observe on Earth. 

Detailed information for the properties of the binary 
population may be obtained through population synthesis. 
We use the binary population synthesis code SeBa to sim- 
ulate the characteristics of the binary population in the 
Galaxy (Portegies Zwart & Verbunt 1996; Portegies Zwart 
& Yungelson 1998). The characteristics of the extragalactic 
population are derived from extrapolating these results to 
the local Universe. 

The paper is organized as follows: in Section 2 we de- 
scribe the population synthesis calculations. Section 3 deals 
with the energy spectrum of a single source followed, in Sec- 
tion 4, by the derivation of the extragalactic backgrounds 
for the different binary populations. In Sections 3 and 4 we 
also give details on the adopted astrophysical and cosmo- 
logical properties of the systems. In Section 5, we compute 
the spectral strain amplitude produced by each cosmologi- 
cal population and investigate its detectability with LISA. 
Finally, in Section 6 we summarize our main results and 
compare them with other previously estimated astrophysi- 
cal background signals. 



2 POPULATION SYNTHESIS MODEL 

In order to characterize the main properties of different com- 
pact binary systems, we use the binary population synthesis 
program SeBa of Portegies Zwart & Yungelson (1998). De- 
tails of the code can be found in (Portegies Zwart & Yun- 
gelson 1998). Here, we simply recall the main assumptions 
of their adopted model B, which satisfactorily reproduces 
the properties of observed high-mass binary pulsars (with 
neutron star companions). 

The following initial conditions were assumed: the mass 
of the primary star mi is determined using the mass function 
described by Scalo (1986) between 0.1 and 100 M@. For a 
given mi, the mass of the secondary star mi is randomly 
selected from a uniform distribution between a minimum of 
0.1 M© and the mass of the primary star. The semi-major 
axis distribution is taken flat in log a (Kraicheva et al. 1978) 
ranging from Roche-lobe contact up to 10 6 R© (Abt & Levy 
1978; Duquennoy & Mayor 1991). The initial eccentricity 
distribution is independent of the other orbital parameters, 
and is S(e) = 2e. 

Neutron stars receive a velocity kick upon birth. Follow- 
ing Hartman et al. (1997), model B assumes the distribution 
for isotropic kick velocities (Paczyhski 1990), 



with u = v/a and a = 600 kms -1 . 

The birthrate of the various compact binaries is nor- 
malized to the Type Il+Ib/c supernova rate (see Portegies 
Zwart & Verbunt 1996). The supernova rate of 0.01 per year 
was assumed to be constant over the lifetime of the galactic 
disc (~ 10 Gyr). 

When computing the birth and merger-rates we account 
for the time-delay between the formation of the progeni- 
tor system and that of the corresponding degenerate bi- 
nary, t 3 . Its value is set by the time it takes for the least 
massive between the two companion stars to evolve on the 
main sequence. For (bh, bh), (ns, ns) and (bh, ns) systems 
T a 50 Myr and it is negligible compared to the assumed 
lifetime of the galactic disc. Conversely, (wd, wd), (ns, wd) 
and (bh, wd) binaries follow a slower evolutionary clock and 
T a can be considerably larger. The cumulative probability 
distribution, P(< T s ), predicted by the population synthe- 
sis code is shown in Fig. 1. For these systems r s can be as 
large as 10 Gyr although all systems are predicted to have 
r s «S 10 Gyr. 

After the degenerate binary has formed, its further evo- 
lution is determined by the time it takes to radiate away its 
orbital energy in gravitational waves. The time between the 
formation of the degenerate system and its final coalescence, 
r m , depends on the orbital configuration and on the mass of 
the two companion stars. The predicted cumulative proba- 
bility distribution is shown in Fig. 2 for the (wd, wd), (ns, 
ns) and (ns, wd) samples. We see from the figure that there 
is a significant fraction of systems which does not merge in 
10 Gyr. For (bh, bh) binaries and mixed systems with one 
black hole companion the population synthesis code predicts 
very long merger times. In particular, all (bh, bh) systems 
appear to have T m greater than 15 Gyr. The reason for these 
large merger times is that binaries with a black hole com- 
panion are characterized by very large initial orbital sep- 
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Figure 1. The cumulative probability distribution for the time 
delay r B (in Myr) between the formation of the progenitor sys- 
tem and the formation of the corresponding degenerate binary 
obtained for the (bh, wd), (ns, wd) and (wd, wd) samples. 
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Figure 2. The cumulative probability distribution for the merger 
time r m (in Myr) is shown for the (ns, ns), (ns, wd) and (wd, wd) 
samples. 
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Figure 3. The probability distribution for the orbital parameters 
of (bh, bh) systems (left panel) is compared to that obtained for 
the (ns, ns) (right panel) population. 



Table 1. The galactic birthrates, Rx,gal, an d merger rates, 
gal' orj tained for each compact binary type X using model 
B of (Portegics Zwart & Yungelson 1998), see text. The rates are 
normalized to the core-collapse supernova rate of 0.01 yr _1 and 
100% binarity. Merger rates are computed after 10 Gyr of the 
evolution of the Galaxy with a constant supernova rate. 



Binary Type X R x ,gaiyr 



arations (see e.g. Fig. 3). In fact, bh progenitors are very 
massive stars and have a very strong stellar wind. For this 
reason they do not easily reach Roche-lobe over-flow and 
rarely experience a phase of mass transfer, which is required 
to reduce the orbital separation of the stars. Unfortunately, 
the evolution (especially the amount of mass loss in the stel- 
lar winds) of high mass stars is rather uncertain (Langer el 
al. 1994). The result that we obtain at least indicates that 
it will be very rare to observe any of these bh mergers. Re- 
cently Portegies Zwart & McMillan (1999) however identi- 
fied a new channel for producing black hole binaries which 
are eligible to mergers on a relatively short time scale. 

In Table [l], we summarize the results for all binary types 
that we have investigated. 
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3 INSPIRAL ENERGY SPECTRUM OF 
SINGLE SOURCES 

Assuming that the orbit of the binary system has already 
been circularized by gravitational radiation reaction, the in- 
spiral spectrum dE/du emitted by a single source can be ob- 
tained using the quadrupole approximation (Misner, Thorne 
& Wheeler 1995). The resulting expression, in geometrical 
units (G=c=l), can be written as, 

dE 7T M 5/3 

du 3 (tti/)V3 > 
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where M = fj, 3 ^ 5 M 2 ^ 5 is the so-called chirp mass, M = 
mi + 7712 stands for the total mass and /j = mi 777,2 /M is the 
reduced mass. 

The frequency v at which gravitational waves are emit- 
ted is twice the orbital frequency and depends on the time 
left to the merger event through, 



(■KU) 8/3 = (lTV max ) 



-8/3 + ^^5/3 {tc _ t) 




(3) 



where t c is the time of the final coalescence and we termi- 
nate the inspiral spectrum at a frequency u max ■ When Post- 
Newtonian expansion terms are included, it is customary to 
consider the inspiral spectrum as a good approximation all 
the way up to ulsco, i-e. the frequency of the quadrupole 
waves emitted at the last stable circular orbit (LSCO) (see 
e.g. Flanagan & Hughes 1998). However, for the purposes 
of our study, we neglect post-Newtonian terms and we set 
the value of v max to correspond to the quadrupole frequency 
emitted when the orbital separation is roughly 3 times the 
separation at the LSCO. 

The value of the orbital separation at the LSCO de- 
pends on the mass ratio of the two stellar components and 
varies between 5M — 6M. The lower limit is obtained in the 
test particle approximation (mi 3> 012), whereas the up- 
per value corresponds to the equal-mass case (mi ~ 7712) 
(see Kidder, Will & Wiseman 1993). If we consider the 
equal-mass limit, which is more conservative for constrain- 
ing the maximum frequency, vlsco ~ 0.022/M and v m ax ~ 
0.19 ulsco- 

For (wd, wd) binaries and binaries with one wd com- 
panion, the maximum frequency, i.e. the minimum distance 
between the two stellar components, is set in order to cut- 
off the Roche-lobe contact stage. In fact, the mass transfer 
from one component to its companion transforms the origi- 
nal detached binary into a semi-detached binary. This pro- 
cess can be accompanied by the loss of angular momentum 
with mass loss from the system and the above description 
cannot be applied. Thus, for closed white dwarf binaries we 
require that the mimimum orbital separation is given by 
r w d{mi) + r wd {m2) where 



r wd {m) = 0.012i? G 



1.44M, 



-2/3 



1.44M, 



2/3 



(4) 



is the approximate formula for the radius of a white dwarf 
from Nauenberg (1972) (see also Portegies Zwart & Verbunt 
1996). 

Consider now sources at cosmological distances. The 
locally measured average energy flux emitted by a source at 
redshift 2 is, 



. dQ dE 



(l + z) 2 c 3 dE e 



4nd L (z) 2 G dv £ 



K1 + *)] (5) 



where dh{z) is the luminosity distance to the source, v = 
v e (1 + z)^ 1 is the observed frequency and the factor c 3 /G 
is needed in order to change from geometrical to physical 
units. Thus, the emission spectrum can be written as, 



dE e M 5/3 
— H l+z)] = - ¥ - (r - ¥r3 



(6) 



where v is the observed frequency emitted by a system at 
time t(z) 



(™r 8 / 3 = (^ max y^(i + zf' 3 + ±p m 5/3 (7) 

[t(z f )+T m -t(z)](l+zf 



and we have written the time of the final coalescence t c = 
t(z c ) in terms of the time of formation t(zf) and of the 
merger-time r m = t(z c ) — t(zf). 



4 EXTRAGALACTIC BACKGROUNDS FROM 
DIFFERENT BINARY POPULATIONS 

Our main purpose is to estimate the stochastic background 
signal generated by different populations of compact binary 
systems at extragalactic distances. 

These gravitational sources have been forming since the 
onset of galaxy formation in the Universe and for each binary 
type X [(ns, ns); (wd, wd); (bh, bh); (ns, wd); (bh, wd); (bh, 
ns)] we should think of a large ensemble of unresolved and 
uncorrelated elements, each characterized by its masses mi 
and rri2 (or M and fi), by its redshift and by its time-delays 
t s and r m [see eqns Q and (Q)]. 

Thus, in order to consider all contributions from dif- 
ferent elements of the ensemble X, we must integrate the 
single source emission spectrum over the distribution func- 
tions for the masses M and fi, for the time-delays and for 
the redshifts. 

The distribution functions for r s and r m depend on the 
binary type X and have been derived from the population 
synthesis code discussed in the previous section. The distri- 
bution function for M can be similarly estimated. 

However, t s , r m and M are not independent random 
variables. In particular, r m and M are correlated because 
M defines the rate of orbital decay, once the degenerate 
system has formed. Thus, for each binary population X , we 
consider the joint probability distribution for r s , r m and A4, 



Px{r s ,r m ,M) dr m dr 3 dM. 



(8) 



Conversely, the redshift distribution function, i.e. the 
evolution of the formation rate for each binary type X, can 
be deduced from the observed cosmic star formation history 
out to z ~ 5. 

In the following subsections, we illustrate the procedure 
we have followed to derive the birth and merger-rates for 
all binary populations and to compute the spectra of the 
corresponding stochastic gravitational backgrounds. 

4.1 Cosmic star formation rate 

In the last few years, the extraordinary advances attained in 
observational cosmology have led to the possibility of identi- 
fying actively star forming galaxies at increasing cosmologi- 
cal look-back times (see, for a thorough review, Ellis 1997). 
Using the rest-frame UV-optical luminosity as an indicator 
of the star formation activity and integrating on the overall 
galaxy population, the data obtained with the Hubble Space 
Telscope (HST Madau et at 1996, Connolly et al. 1997, 
Madau et al. 1998a) Keck and other large telescopes (Stei- 
del et al. 1996, 1999) together with the completion of several 
large redshift surveys (Lilly et al. 1996, Gallego et al. 1995, 



© 2000 RAS, MNRAS 000, §4l| 



Gravitational waves from cosmological compact binaries 5 




Figure 4. The Log of the star formation rate density in units of 
Mqjt - ^^Mpc - 3 as a function of redshift for a cosmological back- 
ground model with fijvf = 1, £7a = 0, ffrj = 50kms — 1 Mpc — 1 and 
a Scalo (1986) IMF. Left: The data points correspond to UV, 
Ha and IR observations of field galaxies. (Filled dots) UV ob- 
servations of Treyer et al. (1998), Lilly et al. (1996), Connolly et 
al. (1997), HDF Madau et al. (1996), Steidel et al. (1996), (1999); 
(asterisks) Ho observations of Gallego et al. 1995, Gronwall 
1999, Tresse & Maddox 1998, Glazebrook et al. 1998, Yan et 
al. 1999); (triangles): ISO IR observations (Flores et al. 1999) 
and the lower limit of SCUBA data (Hughes et al. 1998). (Right): 
The dust-corrected SFR density as derived from UV data in the 
two models most favoured by observations predicted by Calzctti 
& Heckman (asterisks) and by Pei, Fall & Hauser filled dots 
(see text). 



Treyer et al. 1998) have enabled, for the first time, to de- 
rive coherent models for the star formation rate evolution 
throughout the Universe. 

A collection of some data obtained at different redshifts is 
shown in the left panel of Figure ^ for a flat cosmological 
background model with SI a = 0, h = 0.5 and a Scalo (1986) 
IMF with masses in the range 0.1 — IOOA/q. Although the 
strong luminosity evolution observed between redshift and 
1-2 is believed to be quite firmly established, the behaviour 
of the star formation rate at high redshift is still relatively 
uncertain. In particular, the decline of the star formation 
rate density implied by the < z > ~ 4 point of the Hub- 
ble Deep Field (HDF, see Fig. g) is now contradicted by 
the star formation rate density derived from a new ground- 
based sample of Lyman break galaxies with (z) = 4.13 (Stei- 
del et al. 1999) which, instead, seems to indicate that the 
star formation rate density remains substantially constant at 
z > 1 — 2. It has been suggested that this discrepancy might 
be caused by problems of sample variance in the HDF point 



at < z > = 4 (Steidel et al. 1999). 

Because dust extinction can lead to an underestimate of 
the real UV-optical emission and, ultimately, of the real star 
formation activity, the data shown in the left panel of Fig. [| 
need to be corrected upwards according to specific models 
for the evolution of dust opacity with redshift. In the right 
panel of Fig. ^, the data have been dust-corrected according 
to factors obtained by Calzetti & Heckman (1999) and by 
Pei, Fall & Hauser (1999). Using different approaches, these 
authors have recently investigated the cosmic histories of 
stars, gas, heavy elements and dust in galaxies using as in- 
puts the available data from quasar absorption-line surveys, 
optical and UV imaging of field galaxies, redshift surveys 
and the COBE DIRBE and FIRAS measurements of the 
cosmic IR background radiation. The solutions they obtain 
appear to reproduce remarkably well a variety of observa- 
tions that were not used as inputs, among which the SFR 
at various redshifts from Ha, mid-IR and submm observa- 
tions and the mean abundance of heavy elements at various 
epochs from surveys of damped Lyman-a systems. 

As we can see from the right panel of Fig. ^, spectro- 
scopic and photometric surveys in different wavebands point 
to a consistent picture of the low-to-intermediate redshift 
evolution: the SFR density rises rapidly as we go from the 
local value to a redshift between ~ 1 — 2 and remains roughly 
flat between redshifts ~ 2 — 3. At higher redshifts, two dif- 
ferent evolutionary tracks seem to be consistent with the 
data: the SFR density might remain substantially constant 
at z ^ 2 (Calzetti & Heckman 1999) or it might decrease 
again out to a redshift of ~ 4 (Pei, Fall & Hauser 1999). 
Hereafter, we always indicate the former model as 'mono- 
lithic scenario' and the latter as 'hierarchical scenario' al- 
though this choice is only ment to be illustrative. In fact, 
preliminary considerations have pointed out that a constant 
SFR activity at high redshifts might not be unexpected in hi- 
ererachical structure formation models (Steidel et al. 1999). 

Thus, we have updated the star formation rate model 
that we have considered in previous analyses (FMSI, FM- 
SII), even though the gravitational wave backgrounds are 
more contributed by low-to-intermediate redshift sources 
than by distant ones. In addition, if a larger dust correc- 
tion factor should be applied at intermediate redshifts, this 
would result in a similar amplification of the gravitational 
background spectra. 



4.2 Birth and merger rate evolution 

Following the method we have previously proposed (FMSI, 
FMSII), for each binary type X the birth and merger-rate 
evolution could be computed from the observed star forma- 
tion rate evolution. However, this procedure proves to be 
unsatisfactory because it fails to provide a fully consistent 
normalization. Its main weakness is that, even if we assume 
100% of binarity, i.e. that all stars are in binary systems, 
the star formation histories that we have described above 
are not corrected for the presence of secondary stars. For 
the mass distributions that we have considered, secondary 
stars are expected to give a significant contribution to the 
observed UV luminosity as they account for ~ 1/3 of the 
fraction of mass in stars more massive than 8 Mq . 

In order to circumvent the necessity of extrapolating 
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Figure 5. The rest-frame frequency of core-collapse SNae vs red- 
shift predicted by the monolithic and hierarchical models. The 
predictions are consistent with the observed value for the present- 
day galaxy population (see text). 



the UV luminosity indication of massive star formation to 
the full range of stellar masses predicted by the model, we 
could directly normalize to the rate of core-collapse super- 
novae. This is consistent with the adopted normalization for 
galactic rates. 

The core-collapse supernova rate can be directly derived 
from the observed UV luminosity at each redshift, as stars 
which dominate the UV emission from a galaxy are the same 
stars which, at the end of the nuclear burning, explode as 
Type II+Ib/c SNae. Moreover, the supernova rate is ob- 
served independently of the SFR. Therefore it can be used 
as an alternative normalization. 

The rates of core-collapse supernovaepredicted by the 
models shown in Fig. ^] are shown in Fig. ^ assuming a fiat 
cosmological background model with zero cosmological con- 
stant and h = 0.5. 

In the same figure, we have plotted the available ob- 
servations for the core-collapse supernova frequency in the 
local Universe (Cappellaro el al. 1997, Tamman el al. 1994, 
Evans el al. 1989, see also Madau el al. 1998b). 

The binary birthrate per entry per year and comoving 
volume fj(z) can be related to the core-collapse supernova 
rate ny SNaeII+ j b / ( J i (z) shown in Fig. | in the following way, 



fi(z) 



n[SNaeII + Ib/c} ( z ) 
N[SNaeII+Ib/c] 



(9) 



where NysNaeii+ib/c] is the total number of core-collapse 
supernovae that we find in the simulation. 

In order to estimate, from r)(z), the birth and merger- 
rate evolution of a degenerate binary population X, we need 
to multiply eq. ^ by the number of type X systems in 
the simulated samples, Nx, and we also need to properly 
account for both t s and r m . 

We shall assume that the redshift at the onset of galaxy 
formation in the Universe is zf — 5 and that a zero-age main 
sequence binary forms at a redshift z s . After a time interval 
t 3 , the system has evolved into a degenerate binary. Then, 
the redshift of formation of the degenerate binary system, 
Zf, is defined as t(zf) = t(z 3 ) + t 3 . The system then evolves 
according to gravitational wave reaction until, after a time 



interval r m , it finally merges. Thus, the redshift at which 
coalescence occurs, z c , is given by t(z c ) — t[z/) +r m . 

This simple picture implies that the number of X sys- 
tems formed per unit time and comoving volume at redshift 

zf is 



iix(zf) 



dr m dM 



t{z f )~t(z F ) 



n[SNaeII+Ib/c] (Zs) 
(1+Z S ) 



dr a fx 



Px(T s ,T m ,M) 



(10) 



where fx = N x /N[ S N ae ll+lb/c] and z s is defined by t(z a ) 
t(z f ) - Ts. 

If we write, 



px(i 



where r. 



)5{t, 



m i m 



,i)5[M-Mi){ll) 



s,2, ' m,i 



and Mi indicate the time delays and the 
chirp mass for the i th element of the ensemble X, the 
birthrate reads, 



h x (zf) 



1 



N[SN ael I-\-Ib/ c] 

e[t(z f ) -t{ ZF ) - Ts ,i] 



n[SNaeII + Ib/c]{Zs) 



{1 + Zs 



where Q(x) is the step-function. 

Similarly, the number of X systems per unit time and 
comoving volume which merge at redshift z c is, 



t(z c )-t(z F ) 



dr n 



t{z e )- 



-t(z F ) 



10 JO 

ri[SNaeII + Ib/c] ( z s) 



(1+Zs) 



dr s fx 
dM M) 



(13) 



where z s is defined by t(z a ) = t(z c ) — r m — t s . If we apply 
eq. (p. l|) , we can write the merger-rate in a form similar to 
eq. (112), i.e. 



N[SNaeII + Ib/c] 
Q[t(z c ) - t{ ZF ) 



E 



n[SNaelI + Ib/c] ( z s) 
(1 + Zs) 



(14) 



Using this procedure, we compute the birth and merger- 
rates for all the synthetic binary populations. The results are 
presented in Figs. 6, 7 and 8. 

Due to their relatively small t s compared to the cos- 
mic time, the birthrates of (bh, bh), (ns, ns) and (bh, ns) 
systems closely trace the UV-luminosity evolution, although 
with different amplitudes. Our simulation suggests that (bh, 
bh) systems are more numerous than (ns, ns) or (bh, ns) (see 
Fig. 6). 

Conversely, Fig. 7 shows that the birthrates of (wd, 
wd), (ns, wd) and (bh, wd) systems misrepresent the origi- 
nal UV-luminosity evolution as a consequence of their large 
t s . The largest is the characteristic time-delay t s , the more 
the maximum is shifted versus lower redshifts because the 
intense star formation activity observed at z ^ 2, especially 
for monolithic scenarios, boosts the formation of degenerate 
systems at z ^ 2. For hierarchical scenarios, if the redshift 
at which significant star formation begins to occur is z F ~ 5, 
the birthrate of degenerate systems at redshifts ^ 4 is al- 
most negligible. 
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Figure 6. The formation rate of (bh, bh) (asterisks), (ns, ns) 
(dots) and (bh, ns) (triangles) binaries as a function of z is 
shown for hierarchical (upper panel) and monolithic (lower 
panel) scenarios for a flat cosmological background. 
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Figure 7. The formation rate of (wd, wd) (asterisks), (ns, 
wd) (dots) and (bh, wd) (triangles) binaries as a function of 
2 is shown for hierarchical (upper panel) and monolithic (lower 
panel) scenarios for a flat cosmological background. 



Finally, in Fig. 8, we have shown the predicted merger- 
rate for (wd, wd), (ns, ns) and (ns, wd) systems. In this 
case, the distortion of the original UV-luminosity evolution 
is even more apparent, particularly for monolithic scenarios. 
The redshift at which the maximum merger-rate occurs as 
well as the high redshift tail reflects the different r m distribu- 
tions of these populations. We have not shown the merger- 
rates for (bh, bh), (bh, wd) and (bh, ns) binaries because, as 
we have discussed in the previous section, these systems are 
predicted to have merger-rates consistent with zero through- 
out the history of the Universe as a consequence of their very 
large initial orbital separations. 



4.3 Stochastic backgrounds 

Having characterized each ensemble X by the distribution 
of chirp mass and time delays, Px(t s , t to , M), and by the 
birthrate density evolution per entry r/(z), we can sum up 
the gravitational signals coming from all the elements of the 
ensemble. The spectrum of the resulting stochastic back- 
ground, for a binary type X and at a given observation fre- 
quency v, is given by the following expression, 



m&w - I dzfJ o s ~(r+^y (15) 

~° dV 

dM dTmPx(T s ,T m , M) -j—; f[v, Z*] 

where zf is the redshift of the onset of star formation in the 
Universe, Zf is the redshift of formation of the degenerate 
binary systems, z s is the redshift of formation of the corre- 
sponding progenitor system defined by t(z a ) = t(zf) — r a , 
f[v, z%\ is given by eq. (5) and z% is the redshift of emission 
that an element of the ensemble must have in order to con- 
tribute to the energy density at the observation frequency 
v. 

It follows from eq. (f?|) that, for a given observation fre- 
quency u, z* is a function of zj, T m , M and u max . In prin- 
ciple, an inspiraling compact binary system emits a contin- 
uous signal from its formation to its final coalescence thus, 
z c ^ z R ^ Zf. However, in eq. ( |l5| ) we do not restrict to 
systems which reach their final coalescence at z c ^ as we 
are interested to any source between z = and z = zf 
emitting gravitational waves during its early inspiral phase. 
Therefore, the signals which contribute to the local energy 
density at observation frequency v might be emitted any- 
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(dots) and (ns, wd) (triangles) binaries as a function of z 
is shown for hierarchical (upper panel) and monolithic (lower 
panel) scenarios for a flat cosmological background. 



where between sup[0,z c ] ^ z* ^ Zf, provided that, 

/ -1-8/3 / ^8/3/i i *\8/3 , 256 , .5/3 , 1c n 
(TTV) ' = (-RVmax) ' (1 + Z e ) 1 H — M ' (16) 

o 

[t( Zf )+T m -t(z:)](i+z:) 8/3 . 

Substituting eq. ( |ll|) in eq. ([li)), we can write the back- 
ground energy density generated by a population X in the 
form, 



dE 



dYidtdu 



dz 



Nx ■/ i 



(1 + 2.) 



(17) 



dV 



@[t{Zf)-t{zF)-Ts,i]-^f\v, Zl) 



where z% satisfies eq. (|l6|). 

The predicted spectral energy densities for the popula- 
tions of degenerate binary types that we have considered are 
plotted in Fig. 9. For each binary type, we show the results 
obtained assuming both monolithic and hierarchical scenar- 
ios for the evolution of the underlying galaxy population. 

The spectral energy densities are characterized by the 
presence of a sharp maximum which, depending on the bi- 
nary population, has an amplitude spanning about two or- 
ders of magnitudes, in the frequency range [10 _5 -10- 4 ] Hz. 
In the following, we refer to this part of the signal as 'pri- 



mary' component. At higher frequencies, a 'secondary' com- 
ponent appears for all but (bh, bh) systems. The frequency 
which marks the transition between primary and secondary 
components as well as their relative amplitudes depend sen- 
sitively on the population. 

The reason why (bh, bh) systems do not show a sec- 
ondary component is that this is entirely contributed by 
sources which merge before z — 0. Conversely, the low- 
frequency part of the spectrum is dominated by systems 
with merger-times larger than a Hubble time. These sources 
are observed at very low frequencies because the value of 
the minimum frequency (which is emitted at formation, zf) 
is set by the amplitude of the merger-time [see eq. (16)]. 
The larger is the merger-time, the smaller the minimum 
frequency at which the in-spiral waves are emitted. More- 
over, eq. (6) shows that the flux emitted by each source de- 
creases with frequency. This explains the larger amplitude 
of primary components with respect to secondary ones. For 
systems with merger-times larger than a Hubble time, the 
largest frequency is emitted at z = by binaries which form 
at Zf ~ zf- No contribution from such objects can be ob- 
served above this critical frequency and the primary com- 
ponent falls rapidly to zero. 

The amplitude of secondary components reflects the 
number of systems with moderate merger-times. The maxi- 
mum frequency which might be observed is emitted by sys- 
tems which are very close to their coalescence at z = 0. Since 
Vmax is larger for (ns, ns) than for (wd, wd), the secondary 
component produced by double neutron stars extends up to 
~ 10 2 Hz. 

It is interesting to note is that monolithic scenarios pre- 
dict a maximum amplitude which is a factor ~ 20-25% larger 
than the hierarchical case. This difference is much larger 
than what has been previously obtained for other extra- 
galactic backgrounds (see e.g. FMSI), indicating that the 
energy density produced by extragalactic compact binaries 
is substantially contributed by sources which form at red- 
shifts ^ 1 — 2. It is quite difficult to unveil the origin of this 
effect because of the large number of parameters which de- 
termine the appearance of the final energy density. However, 
a plausible explanation might be that, depending on its spe- 
cific time-delays t s and r m , each system emits the signal at 
redshifts which can be substantially smaller than the forma- 
tion redshift of the corresponding progenitor system. Thus, 
although the background signal is mostly emitted at low- 
to-intermediate redshifts, the sources which produce these 
signals might have been formed at higher redshifts and re- 
flect the state of the Universe at earlier times, when the 
differences among hierarchical and monolithic scenarios are 
more significant. Comparing the different panels of Fig. 9, we 
conclude that the background produced by (bh, bh) binaries 
has the largest amplitude but it is concentrated at frequen- 
cies below ~ 2 x 10~ 5 Hz. At higher frequencies, which are 
more interesting from the point of view of detectability, the 
dominant contribution comes from (wd, wd) systems. This is 
consistent with what has already been found for the galactic 
populations (Hils, Bender & Webbink 1990). 

^From the background spectrum it is possible to com- 
pute the closure density Q. gw h 2 and the spectral strain am- 
plitude of the signal Sh, 



dE 



7tc 3 v 2 dSdtdv 



(18) 
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Figure 9. The spectral energy density of the gravitational background produced by various extragalactic populations of degenerate 
binaries in monolithic and hierarchical scenarios assuming a flat cosmological background with zero cosmological constant. 



= i^ir^rM- (19) 

c 3 p cr dtdSdu 

The results are shown in Figs. 10 and 11 for all binary types 
within monolithic and hierarchical scenarios. 

The strain amplitude of the backgrounds has a max- 
imum amplitude between ~ 10 -18 Hz -1 / 2 and ~ 5 x 
10 -1 ' Hz" 1 / 2 at frequencies in the interval [~ 5 x 10~ 6 — 
5 x 1CT 5 ] Hz. The function S h is more sensitive to the low 
frequency part of the energy density. Therefore, its shape re- 
flects mainly the primary components of the corresponding 
energy density. In all but the (bh, bh) population, it is evi- 
dent the presence of a tail at frequencies above the maximum 
which is the secondary component of the energy density: in 
the next section we compare this part of the background 
signal with the LISA sensitivity to assess the possibility of 
a detection. Still, it is clear that the prominent part of the 
background signals produced by extragalactic populations 



of degenerate binaries could be observed with a detector 
sensitive to smaller frequencies than LISA. 

Conversely, Q gw h 2 is mostly dominated by secondary 
components. We can compare the predictions for (bh, bh), 
(wd, wd) and (ns, ns) systems. Contrary to what has been 
found for the spectral energy density or for the strain am- 
plitude of the signal, the largest Q gw h 2 is produced by (ns, 
ns), as a consequence of the high amplitude of the secondary 
component. In particular, no significant contribution from 
the primary component appears. For (wd, wd), instead, the 
contribution of the primary component is relevant, although 
its amplitude is roughly half that of the secondary compo- 
nent. Finally, for (bh, bh) no secondary component is pro- 
duced and thus the amplitude of the closure density is very 
low and at very low frequencies. Mixed binary types have 
different properties, depending on the relative importance of 
the above effect. For instance, (bh, wd) produce a secondary 
component but the amplitude is so small to be comparable 
with that of the primary. 
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Figure 10. The strain amplitude of the gravitational background produced by various extragalactic populations of degenerate binaries 
in monolithic and hierarchical scenarios assuming a flat cosmological background model with zero cosmological constant. 



We stress that the value of v max is quite uncertain as it 
defines the boundary between the early inspiral phase and 
the highly non-linear merger. Clearly, the more we get closer 
to this boundary, the less accurate is the Newtonian descrip- 
tion of the orbit as post-Newtonian terms start to be rele- 
vant. Therefore, we believe that the most reliable part of the 
binary background signal is the low frequency part, i.e. the 
part which mostly contributes to the strain amplitude Sh- 



5 CONFUSION NOISE LEVEL AND 
DETECTABILITY BY LISA 

To have some confidence in the detection of a stochastic 
gravitational background with LISA it is necessary to have 
a sufficiently large SNR. The standard choice made by the 
LISA collaboration is SNR = 5 which, in turn, yields a mini- 
mum detectable amplitude of a stochastic signal of (see Ben- 
der 1998 and references therein), 



h 2 VL gw [v = lmHz] ~ KT 12 . (20) 

This value already accounts for the angle between the 
arms (60°) and the effect of LISA motion. It shows the re- 
markable sensitivity that would be reached in the search for 
stochastic signals at low frequencies. Table |^ shows that the 
backgrounds generated by (wd, wd) and (ns, ns) extragalac- 
tic binary populations exceed this minimum value and LISA 
might be able to detect these signals. 

We plot in Fig. ^ the predicted sensitivity of LISA to 
a stochastic background after 1 year of observation (Bender 
1998). On the vertical axis it is shown forms, denned as, 

/ ir ms = [2i/S, l ( I z)] 1/2 (^) V2 (21) 

where S„(v) is the predicted spectral noise density and the 
factor (Aiz/j/) 1 / 2 is introduced to account for the frequency 
resolution Au = 1/T attained after a total observation time 
T. 
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Figure 11. The function £l gw h 2 of the gravitational background produced by various extragalactic populations of degenerate binaries 
in monolithic and hierarchical scenarios assuming a flat cosmological background model with zero cosmological constant. 



v =1 mHz (wd, wd) 
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Table 2. The values of the closure density at 1 mHz obtained 
for (wd, wd) and (ns, ns) extragalactic binary populations inves- 
tigated. All these values are larger than the minimum detectable 
value of Q gw /i 2 (lmHz) predicted by the LISA team for a SNR = 5 



and after 1 year of observation. 



On the same plot we show the equivalent forms levels pre- 
dicted for different extragalactic binary populations and for 
the galactic population of close white dwarfs binary consid- 
ered by Bender & Hils (1997). 

We see that the extragalactic backgrounds might be ob- 
servable at frequencies between ~ 1 and ~ 10 mHz. 



These background signals represent additional noise 
components to the LISA sensitivity curve when searching 
for signals from individual sources. 

In particular, backgrounds from unresolved astrophysi- 
cal sources represent a confusion limited noise. In fact, unless 
the signal emitted by an individual source has a much higher 
amplitude, the background signal prevents the individual 
source to be resolved. Clearly, the magnitude of this effect 
depends on the frequency resolution of the instrument, i.e. 
on the observation time. The forms noise levels produced by 
extragalactic compact binaries shown in Fig. [l^ have been 
computed assuming T = 1 yr. For the same total observa- 
tion time we show, in Fig. |l^, the number of extragalactic 
(wd, wd) and (ns, ns) observed in each frequency resolution 
bin. At frequencies were these backgrounds might be rele- 
vant (between 1 and 10 mHz), the number of sources per 
bin is 3> 1, representing a relevant confusion limited noise 
component. The critical frequency above which the number 
of sources per bin is lower than 1 occurs at ~ 0.1 Hz for (wd, 
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Figure 12. The sensitivity of LISA to a stochastic background 
of gravitational waves after one year of observation. The extra- 
galactic backgrounds from (wd, wd), (ns, ns), (ns, wd) and (bh, 
ns) systems might be observable at frequencies between ~ 1 and 
~ 10 mHz. 
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Figure 13. The number of extragalactic (wd, wd) and (ns, ns) 
binaries per resolution bin after a total observation time of lyr. 

wd) and outside LISA sensitivity window for (ns, ns). How- 
ever, at these frequencies the dominant noise component is 
the instrumental noise. 



6 CONCLUSIONS 

In this paper we have obtained estimates for the stochastic 
background of gravitational waves emitted by cosmological 
populations of compact binary systems during their early- 
inspiral phase. 

Since we have restricted our investigation to frequen- 
cies well below the frequency emitted when each system ap- 
proaches its last stable circular orbit, we have characterized 
the single source emission using the quadrupole approxima- 
tion. 

Our main motivation was to develop a simple method 
to estimate the gravitational signal produced by populations 
of binary systems at extragalactic distances. This method 
relies on three main pieces of information: 

(i) the theoretical description of gravitational waveforms 



to characterize the single source contribution to the overall 
background 

(ii) the predictions of binary population synthesis codes 
to characterize the distribution of astrophysical parame- 
ters (masses of the stellar components, orbital parameters, 
merger times etc.) among each ensemble of binary systems 

(iii) a model for the evolution of the cosmic star formation 
history derived from a collection of observations out to z ~ 5 
to infer the evolution of the birth and merger rates for each 
binary population throughout the Universe. 

As we have considered only the early-inspiral phase of the 
binary evolution, our predictions for the resulting gravi- 
tational signals are restricted to the low frequency band 
1CP 5 — 1 Hz. The stochastic background signals produced 
by (wd, wd) and (ns, ns) might be observable with LISA 
and add as confusion limited noise components to the LISA 
instrumental noise and to the signal produced by binaries 
within our own Galaxy. The extragalactic contributions are 
dominant at frequencies in the range 1 — 10 mHz and limit 
the performances expected for LISA in the same range, 
where the previously estimated sensitivity curve was attain- 
ing its mimimum. 

We plan to extend further this preliminary study and to 
consider more realistic waveforms so as to enter a frequency 
region interesting for ground-based experiments. 

Finally, in Fig. [14 we show the spectral densities of the 
extragalactic backgrounds that have been investigated so 
far. The high frequency band appears to be dominated by 
the stochastic signal from a population of rapidly rotating 
neutron stars via the r-mode instability (see FMSII). For 
comparison, we have shown the overall signal emitted during 
the core-collapse of massive stars to black holes (see FMSI). 
In this case, the amplitude and frequency range depend sen- 
sitively on the fraction of progenitor star which participates 
to the collapse. The signal indicated with BH corresponds to 
the conservative assumption that the core mass is ~ 10% of 
the progenitor's (see FMSI). Recent numerical simulations 
of core-collapse supernova explosions (Fryer 1999) appear to 
indicate that for progenitor masses > 40Mq no supernova 
explosion occurs and the star directly collapses to form a 
black hole. The final mass of this core depends strongly on 
the relevance of mass loss caused by stellar winds (Fryer & 
Kalogera 2000). If massive black holes are formed the re- 
sulting background would have a larger amplitude and the 
relevant signal would be shifted towards lower frequencies, 
more interesting for ground-based interferometers (Schnei- 
der, Ferrari & Matarrese 1999). 

In the low frequency band, we have plotted only the 
backgrounds produced by (bh, bh), (ns, ns) and (wd, wd) 
binaries because their signals largely overwhelm those from 
other degenerate binary types. 

We find that both in the low and in the high fre- 
quency band, extragalactic populations generate a signal 
which is comparable to and, in some cases, larger than 
the backgrounds produced by populations of sources within 
our Galaxy (Giazotto, Bonazzola & Gourgoulhon 1997; Gi- 
ampieri 1997; Postnov 1997; Hils, Bender & Webbink 1990; 
Bender & Hils 1997; Postnov & Prokhorov 1998; Nelemans, 
Portegies Zwart & Verbunt 1999). It is important to stress 
that even if future investigations reveal that the amplitude 
of galactic backgrounds might be higher than presently con- 
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Figure 14. The predicted strain amplitude of the stochastic 
backgrounds produced by extragalactic populations of gravita- 
tional sources. In the high frequency band, we show the esti- 
mates for the background produced by rotating neutron stars via 
r-mode instability, and two possible signals emitted by popula- 
tions of massive stars collapsing to black holes (see text). In the 
low frequency band, we plot the background predicted for three 
different populations of binary systems. 



ceived, their signal could still be discriminated from that 
generated by sources at extragalactic distance. In fact, the 
signal produced within the Galaxy shows a characteristic 
amplitude modulation when the antenna changes its orien- 
tation with respect to fixed stars (Giazotto, Bonazzola & 
Gourgoulhon 1997; Giampieri 1997). 

The same conclusions can be drawn when the extra- 
galactic backgrounds are compared to the stochastic relic 
gravitational signals predicted by some classical early Uni- 
verse scenarios. The relic gravitational backgrounds suffer 
of the many uncertainties which characterize our present 
knowledge of the early Universe. According to the presently 
conceived typical spectra, we find that their detectability 
might be severely limited by the amplitude of the more re- 
cent astrophysical backgrounds, especially in the high fre- 
quency band. 
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